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' We study a model of the quantum critical point of cuprates associated with the "circulating 

' current" order parameter proposed by Varma. An effective action of the order parameter in the 

quantum disordered phase is derived using functional integral method, and the physical properties of 
f~| ' the normal state are studied based on the action. The results derived within the ladder approxima- 

I tion indicate that the system is like Fermi liquid near the quantum critical point and in disordered 

regime up to minor corrections. This implies that the suggested marginal Fermi liquid behavior 
' induced by the circulating current fluctuations will come in from beyond the ladder diagrams. 

PACS numbers: 71.10.-w, 71.10.Hf, 74.20.Mn 
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■ I. INTRODUCTION 

c/2 ' In spite of the intensive studies of the last 15 years the microscopic origin of the high temperature superconduc- 
tors (HTSG^ still remains elusive. However, the physical properties in each regime of phase diagram are now well 
^ ' established.u The normal state properties of the cuprate supercondcutors are highly anomalous. They deviate sub- 
stantially from the Fermi liquid behavior and are believed to be of a non-Fermi liquid. Especially, the anomalous 
' ■ normal state properties near optimal doping concentration can be characterized by the phcnomcnological marginal 
^ I Fermi liquid (MFL) theory.H The MFL theory rests on a single hypothesis that there exist excitations whose spectral 
Q r function is given by 
O 



> 



ImP(c., q) ^ <; » (1) 



where N{Ep) is the unrenormalized density of states at the Fermi energy, and sgnw — uj/\uj\. The important features 
1 of the spectral function Eq. (|l|) are its independence of the momentum q and its scaling form as a function of uj/T. 
' The oj/T scaling is the hallmark of quantum critical point (QCP), thus the scaling form of the spectral function very 
strongly suggests a proximate QCP. If QCP indeed underlies the MFL behaviour, the momentum independence of 
Eq. (|l|) implies infinitely large dynamical exponent Zd of QCP. 

Most of the normal state properties associated with two-particle correlations near optimal doping have been found 
' consistent with the MFL theory.H Recently, the validity of MFI^was confirmed also for single-particle properties. The 
recent angle- resolved photoemission experiment by Valla et alB found that the single-particle self-energy is given by 
the MFL form. They measured the imaginary part of the electron self-energy of the optimal Bismuth HTSC along a 
nodal direction, and obtained ImS(k, e) ~ max(|e|, T). This observation renewed the interest in the MFL theory, and 
. especially, in the microscopic basis of the MFL and conjectured proximate QCP. np 
' There are several models rooted on the idea of QCP which attempt to explain various aspects of HTSChEI. In 
Q , general, it is necessary to introduce a certain local order parameter to define a QCP. The present experimental data 
O ■ put strong constraint on the possible form of the order parameter. In particular, there is no evidence of broken 
[ translational symmetryp.|ao the order parameter of QCP should be chosen to respect the translational symmetry. The 
• '~j ■ theory by C. M. VarmaoB is based on the order parameter which does respect the translation symmetry but breaks the 
' time reversal symmetry and the four-fold rotation symmetry. The construction of such order parameter necessitates 
^ . the full set of the copper d and the oxygen p^.y orbitals. The explicit form of the (complex) order parameter is 

*(T, = I 5] (sin ^4.p.p. - sin ^^.p.p.), (2) 

p,(T 

where V is coupling constant, and a, a; is lattice constant of unit celfl and doping concentratio|H_| respectively, a is 
spin quantum number. At zero temperature, '^{T = 0,x < Xc) ^ 0. The critical doping Xc ^ O.20a defines the QCP. 

In QCP theory of HTSC by Varma, the anomalous normal state corresponds to the finite temperature region of 
the disordered regime T > Tc > 0, a; ~ Xc- In this regime the critical fluctuations of order parameter determine the 
physical properties. Varma made an extensive study on the general class of models which can give rise to the order 



parameter (||), and gave an heuristic (and essentially non-perturbative)|Xlerivation of the susceptibility q) which 
controls the critical fluctuation of order parameter in the normal stateQ. In this paper, we have chosen a simplified 
model considered by Varmafl and derived an effective action (See Eg. (pi]).) of the order parameter systematically in 
the functional integral formalism. Such effective action can be also regarded as time-dependent Landau- Ginzburg type 
free energy. The lowest order term of the effective action (See Eq.(p7|).) amounts to the summation of the ladder- type 
diagrams. The physical properties deduced from the lowest order term of the effective action are not compatible with 
the experimental data of the normal states. For an example, the electron self-energy derived from the lowest order 
effective action is Fermi liquid-like apart from logarithmic correction. If we accept the heuristic derivations by Varma, 
the inadequacy of the lowest order approximation implies the importance of higher order terms beyond the ladder 
diagrams. 

This paper is organized as follows: In Sec. II, we introduce our model and derive the functional integral representation 
of the effective action. In Sec. Ill, the concrete form of the effective action is obtained up to the 4th order. In Sec. IV, 
we analyze the effective action following Moriya's self-consistent renormalization scheme, and calculate some physical 
quantities. In Sec.V, we conclude this paper with the summary and some concluding remarks. 



II. FORMULATION 

The basic model is the three-band Hubbard modeliflR 



H — Hq + Hi + H2 , 



H2 = 2V c^in)dl^^^d]^^pl^,^^^,p^k'a' +Cy{q)dl^^^dkaplk'-qa'Py^'<^'- (3) 

k,k' ,q,a",cr' 

I is a lattice site index, and Cx,y(k) = cos -s^^ Sx^yCk) = sin !^s^. Hi is the on-site repulsion term which suppresses 
the doubly occupied states. H2 is the nearest neighbor repulsion {V > 0) term between copper d and oxygen Px,Py 
orbitals written in momentum space. The order parameter of circulating current (CC) phase comes from i?2- To 
simplify our problem, we implement the effect of Hi by a substitution tpd — > xtpd, tpp xtppEl. Certainly, this is very 
crude approximation and the effect of Hi shoud be more carefully included in due course. 
It is easy to show, by re-indexing momenta, that H2 can be written as 



H2 


= -y E 

i=l,2,3,4(T(T' 


, 




k' 


' - Sk'y^k', 


02a'a{q) 


= (sk':rpIk'<T 
k' 


' + Sk'y^k', 




k' 


■ ^ Ck'ypj^k', 




k' 


■ + Ck'y^k', 



C?k'+qo-- (4) 

The expectation value of the operator Oi^a = Oi^ is the order parameter of CC phase Eq.(^) up to constant factor. 
Because the critical fluctuation of Oia is most dominant near QCP, the remainder H'^ = -VJ2^=2,3A,a'aOlfi,^ - 
V^^i_^^0\i^^Ocycyii can be treated as perturbations near QCP. We will not take the effect of H2 into account 

explicitly in this work. Note that the operator Oia is odd under the C4 symmetry of square lattice. The reduced 
Hamiltonian which incoporates the physics associated with CC phase is 

-ff = ^ (ed Ridker + 2tdj^^ [Sa;(k)pa;kcr + Sy{\^)py]^„] - U' Sx{'k)Sy{'k)pl^^^Pyka- + H.C. 
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-v^Eoi(q)o.(q), (5) 

cr,q 

where t = xtpdyt' = xtpp. We could have adopted a basis where kinetic term of the reduced Hamiltonian is diagonal- 
ized. The advantage of our choice of basis is that the couphng between order parameter and fermions is simple and 
that it admits the formal manipulations of functional integral method. Using Hubbard-Stranovich transformation, 
the partition function can be written in the functional integral form. 



,k,k' 



where the six component spinor is defined ^;(j(je,k) — ^dkcnPKkcrjPakcry I, in are the orbital indices. The explicit 
form of kernel matrix is (/i is the chemical potential.) 

/ed-ie- fJ- 2tsx(k) 2tsy{k) \ 

Kq ^ {ie, ie' ■,k,k' ; Im; aa') = See' Saa'Skki i '2tSa:{k) -ie - ^ -4:t'sj;{k)sy{k) \ . (7) 

V 2tsylk) -4fsx{k)sy{k) -ie - fi ) 

The explicit form of the matrix M is (fc = (ie, k), k' — {ie' , k').) 

/ Nk'-ks^ik) ~Nk'-kSyik)\ 

M{k',k;a'a;lm)^6aa' i ivLfc"Sx(k') . (8) 

V-^Lfc'S,(k') / 

The effective action of order parameter can be obtained by integrating out the electrons ^'o-;. {q — {iuj,q).) 

5eff = Vj2Nl{q)NAq) - 2 lnDct(/^o-i - VM), (9) 

aq 

where the factor of 2 of the second term comes from the spin degeneracy. From now on the spin index of N will 
be suppressed. By expanding the determinant of the (time-dependent) Landau- Ginzburg type action can be 
obtained. In principle, the action (|^) contains all the elements for the description of the critical properties of our 
model. The mean-field phase diagram can be determined by minimizing with respect to N,N^. Using relation 
(51ndet A = STihiA = TrA'^M, we get 

Eq.(|lO|) is the mean- field equation written in matrix form. According to the mean- field solution by Varmai, Eq.(|lO|) 
has a non-trivial solution for x < Xc ~ 0.2 and T < Tc{x) with Tc{xc) = 0. The non-trivial mean-field solution is of 
the form N ~ iR, where R is real. It should be noted that only the solution of this type preserves C4 x T, where 
T is the time reversal. Thus, we expect that the susceptibility of the imaginary component of N {l'niN{q)lmN{—q)) 
would diverge as q —> 0,x —> Xc + at zero temperature, reflecting the large critical fluctuation near QCP. The 
susceptibility determines the critical properties near QCP. In the quantum disordered phase, the mean-field equation 
( p^ ) has a trivial solution = 0, and we may expand the action (^) as a power series of N. 



III. EFFECTIVE ACTION 



One can show that the terms with odd powers of A^ vanish in the expansion of the determinant of (||) by using the 
transformation property of A^ under C4 rotation. Expanding the determinant up to the 4-th order we get 

5eff - 2VyN^q)N{q) + y [nA{q)Niq)N^q)+IlB{q) {N {q)N {-q) + N^q)N^ {-q)) 
q 1 

+ y [^MNiqi)N{q2)N^{q3)N^iq4) + n^fe) Niq,)N^{q2) (A(g3)iV(94) + N^{q3)NHq4)) 

(ji, 1=1,2, 3, 4 

+ nE{q^) iNiqi)N{q2)N{q3)Niqi) + NHqi)N\q2)N^q3)NHqi)) \ + OiN^,N^ . . .). (11) 
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The appropriate momentum-energy conservation rules are to be understood for the 4th order terms. The exphcit 
expressions of the polarization functions 11^ are relegated to the appendix. Since the d, Px , Py electrons are the linear 
combinations of the anti-bonding a, the bonding b, and the non-bonding c electron operators, the polarization functions 
Ili's contain both intra-band (aa) and inter-band (ab or ac) contributions. 

For analyses later, it is convenient to decompose the order parameter TV into the real and imaginary components. 
In real space {x = (t, x).), the decomposition is 



N{x) = Ni{x)+iN2{x),N''{x) = Ni{x) -iN2{x). 

In momentum-frequency space, 

Niq) = 7Vi(g) + iN2iq),N\q) = Ni{-q) - iN2{^q). 
In terms of real components Ni , N2 the quadratic part of the effective action (^l]) can be written as 

sis = E [Mq)Ni{q)N,{-q) + ^2{q)N2{q)N2{-q) + tM<l){Ni{q)N2{-q) - N,{-q)N2{q)) 



(12) 



(13) 



(14) 



where 



ir,{q) = 2V + nA{q) + 2IlBiq), 
n2{q)^2V + UA{q)-2IlBiq), 
iT3{q) = 2V + nAiq). 



Using the results Eq.( Bll ) for 11^(9) and IlB{q), 7ri((7)'s {i = 1, 2, 3) can be parametrized as 



TTiiiuj, q) = Si + ai\qf + b^uj'^ + 



(15) 



(16) 



where oi = 03, Fi = F3. The polarization functions 7ri(g) contain both the intra-band and the inter-band contributions. 
In particular, the Landau damping term which exists for \uj\ < Fi|q| solely stems from the intra-band contribution. 
The parametrizations E q.(|l6| ) are valid for |ti;| < fi (See Appendix B.). At this point, we observe that Ni{q)Ni{—q) 
and N2{q)N2{~q) of Eq.(|^are even in q, while {Ni{q)N2{—q) — Ni{—q)N2{q)) is odd in q. Therefore, only the even 



parts of TTi {q) and 7r2 {q) and the odd part of t:^ (q) contribute to S^g . 

•^cff - E [Mq)Ni{q)mi-q)+7T2iq)N2iq)N2i-q) + i7rsiq){N,iq)N2{-q) - TVi (-g) A^2(g)) 



(17) 



T:i{iuj,q) =di+ ai|qp + 

7r2(iw,q) ^62 + a2\qf 
TTaiiuj, q) = -irjuj. 



CO 



ri|q| 

1^1 
r2|q| 



+ blLU^, 
+ b2UJ^, 



The propagators (or the susceptibilities) of order parameters can be read off from Eq.(17jT^). 

■^2{q) 



{Niiq)Nii^q)) = 
{N2iq)N2i-q)) = 
{N,iq)N2i-q)) = 



7i"i((7)7r2((7) - (7r3(g))2' 
Ml) 

-»7^3(g) 
■^i{q)T^2{q) - (T^siq)) 



2 ■ 



Let us define Dij{q) ~ {Ni{q)Nj{—q)) . The explicit form of the Z?22(iw,q) is 



-D22(jt^,q) 



S2 + a2|qp + b20J^ 
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1-1 



r2|q| 5i + ai|q|2-|-6icj2 



Tllql 



(18) 



(19) 



(20) 
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Eq.(pO|) is valid within the ladder approximation (See below.). As discussed in the previous section, the critical 
fluctuation of A'2 is important near QCP, thus the susceptibility 1^22 (*w,q) determines the physical properties near 
QCP. Obviously, 5^ is the parameter which controls the proximity to QCP. The precise location of QCP is determined 
by the condition 



[D(r']-i(zu;^0,q^O) = 0, 



(21) 



where D^2^^\'^^7 q) is the fully renormalized susceptibility which includes all of the higher order corrections. Since Si 



is expected to be large, we can approximate the last term of (^0|) as jj-ijJ^ and absorb it into 62- Then, within ladder 
approximation, the imaginary part of the retarded 022(^,(1) becomes 



IinD§2iio > 0,q) 



e(r2|q 



(^2^ 

r2|q|)^(^2 



«2|qP)^ + (rIR 
-aalqp - 62^^), 



(22) 



where Q{x) is the step function (6(x) = 1, a; > 0; 9(x) = 0, a; < 0.). 

If we neglect the higher order terms like N^^ N^, . . ., the quadratic action Eq.(p7|) amounts to the ladder approxima- 
tion. This can be verified by the direct calculation of susceptibility (O(g)O^(g)) (See Eq.(|4|).) both by the funotional 
integral formalism and the conventional perturbation theory. Then, following Mahan's treatment of excitonsEj, we 
might expect the excition-like poles in the susceptibility. This pole can be identified in Eq. ( pO| ) by analytic continuation 
iuj Lu + i6. QCP in the ladder approximation is equivalent to the vanishing of the exciton-like gap. 

The coefficients of the 4th order terms Ilc{qi),TlD(qi),IlE{qi) depend on the momenta of order parameters. For 
small gi's, Ilc{qi),TlDiQi),^Eiqi) can be expanded in g^. Neglecting the dependences on momenta, the 4th order 
terms of the effective action can be expressed as (in real space x = (r, x).) 



S. 



(4) 



^ Ntix) + I Niix) + A nKx)NI{x) 



(23) 



where 



Ai =4!(nc(0) + 2nc(0) + 20^(0)), A2-4!(nc(0)-2n,,(0) + 20^(0)), A3 = 2!2!(2nc(0) - 120^(0)). (24) 

S'^g introduces the lowest order corrections beyond ladder diagrams. It also renormalizes the electron self-energy at 
one loop order. The importances of still higher order terms will be discussed at Sec.V. 

Near QCP 82 — 0, the naive perturbation theory breaks down because of the infrared divergences as can be 



demonstrated by caculating the correction to free energy with S-p S^'^ + S^g . 
the self-consistent renormalization scheme developed by MoriyaO. 



These divergences can be handled by 



IV. SELF-CONSISTENT RENORMALIZATION 



To cure the infrared divergences we have to determine the critical parameter 82 self-consistently includtog-the effect 
of the 4th order terms Eq.(E3). The self-consistent equation in real space for the renormalized (5^(T) isollj 

5^{T) =62 + ^D22{0;6i,S^{T)) + ^Dn{0;6i,S^{T)), (25) 

where the dependences on the critical parameter are indicated explicitly. We can neglect the renormalization of 61 
since it does not control the critical behaviour. Let 61,62 be the bare value of 61,62 for which 6^{T = 0) = 0. In 
other words, 6^, 62 defines the zero temperature QCP. Then, the temperature dependence of the (5|*'(T) for 6^, 6^ can 
be determined. By subtracting the Eq.(p5|) at T 7^ and T ~ Q ior 61,62, we get 

S^{T) = ^[d22{0;61,6^{T)) ~ D22{0;6l,6i = 0)) + ^[Dn{0;6l,6^{T)) - Dn{Q;6l6^ = 0)). (26) 
Using the relation, 

^'^^°) = J J^l (27) 
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we can reexpress Eq.(E6|) as 



cPq 



duj 
2^ 



A3 

4 



(coth — - l)Imi:)22_5«(-r)(^^,q) 



(28) 



The first term of Eq.(^8|) determines the temperature dependence. The last term of Eq.(^8|) part ) gives the 

subleading contribution proportional to T^, while the first term and the second term give nn((5|*(T) +T^) + const, x 
62 (T). At low temperature, the dominant contribution to the integral comes from the low frequency region, where 
the Landau damping gives the largest contribution to ImD^(ti;, q). Carrying out integrals up to logarithmic accuracy, 
we get 



Tlni 
T 



(29) 



Now we compute some physical properties based on the self-consistent renormalization scheme. 

Specific Heat- The singular behaviour of specific heat near QCP can be obtained by computing free energy with 
renormalized 82 {T). Integrating out order parameter N with the effective action Eq. ([l7|), we get the singular piece 
of the free energy in (self-consistent) ladder approximation. 



(2^)2 



duj uj 

— coth — — Im In 

TT 2r 



7rf(u;,q)^2«(c.,q)-(7r3^(a;,q))^ 



(30) 



The most dominant contribution to the free energy at low temperature comes from the low-w region. In that region, 
we can neglect the terms of polarization functions nf' in the free energy compared to Landau damping term. Then, 
the most singular contribution comes from 7r^(ti;,q). 



smg 

Using the expansion tan^^ y ^ y 
singular part of the specific heat. 



(2^ 



^coth^tan-i fW^^ 
27r 2T \5f{T) 



a2|q|- 



(31) 



j/'^/3, and taking derivative of Ai^sing with respect to temperature we get the 

y3 



T 



(32) 



The above result is not compatible with experimental data in normal state which shows almost T-linear behaviour. 
We have to keep in mind that the result Eq.(^2|) is essentially due to the intra-band contribution to the polarization 
function tt^ . In three dimension, the above result would be 



ji3 



In- 



T 



(33) 



Electron Self-Energy- The coupling between order parameter and the anti-bonding electron (conduction electron) 
can be obtained by substituting Eq. ([A^) into Eq.(||). 



Va 



Ni{k' - fc)/(fc, k') + iN2{k' - k)g{k, k') 



ak, 



where the coupling matrix elements are(notations are defined in Appendix A.) 

2t 



f/g^ 



i^A)k{^A)k' 



EAmsiik) - sUk)) ± EA{h){si{k') - .^(k')) 



(34) 



(35) 



Note that when the momentum transfer k' — k is small the matrix element g{k, k') is proportional to (k — k') • (k + k'), 
therefore, the coupling with the critical order parameter 7V2 is suppressed for the forward scatterings. Considering 
the scattering with N2, the imaginary part of the electron self-energy can be written as 
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X (5(e + w-Cp+q)(^^i3(t^)+^^^^(e + t^)) + 5(e-w-ep+q)(l + ?^s(^^)-?^F(e-w)) , (36) 

where the unimportant factors of the couphngs are omitted. At zero temperature {S2 — 0), the self-energy (|3^ ) 
becomes (assuming e > 0), 



ImS«(e, p) = v'J-^J^^ [(p + |) • cirimDi^iu;, q; = 0)6ie - - ^, 



p+q)- (37) 



Let us first consider the zero temperature case. ImZ3^(ti;, q; = 0) behave differently depending on the relative 
magnitude of iv and r2|q| (see Eq.(^2|).). If w > r2|q|, then the typical values of momentum and frequency are given 
by w ^ g|| ^ q±. q\\/± is the component of q parallel/perpendicular to p. Therefore, the |qp term in the vertex can 
be neglected compared to p • q. Carrying out the integrals, we find the self-energy is proportional to e^. If < r2|q|, 
LV ~ (7||, while qj_ ^ uj^^^. Therefore, qj_ » gy, the p • q term in the vertex can be neglected compared to |qp. 
Carrying out integrals, with the momentum cut-off of the order of the external momenta p, we find the self-energy is 
proportional to —e^lne^^^. The contribution from the region < r|q| is larger than that of from \uj\ > r|q| by a 
logarithmic factor. Thus, the self-energy at zero temperature is 

Iml]^(e,p,r = 0) - -e^lnlel^/^ (38) 

Because the vertex [(p -I- ^) ■ q]^ suppresse the low momentum processes very strongly, the self-energy does not show 
anomalous behaviour except for the logarithmic factor. At finite temperature, the bose factor 71^(0;) dominates the 
integral. Therefore, the frequency integral is effectively cut off by the temperature. Carrying out the integrals, we find 
the result is proportional to —T'^lnT^^^. Combining the results for the zero temperature and the finite temperature 
cases, we can write 

Iml]^(e, p, T) (max(e, T))^ ln[max(e, T)]^/^ (39) 

The self-energy Eq.(^^ is like that of Fermi-liquid apart from the minor logarithmic factor. Eq.(^9|) is also essentially 
determined by intra-band contribution to the polarization function like specific heat. 

The results obtained within self-consistent renormalizaton Eq. (^, |3^) are not compatible with those of MFL theory. 
Clearly, we have to incorporate the higher order corrections beyond ladder approximation. However, the results Eq.(|3^, 
|39| ) might be applicable to overdoped region, namely the deep in the disordered phase, where the critical fluctuation 
is less important, so that the higher order correction is not important. In this region, the temperature dependence of 
62 {T) can be ignored, and the results Eq.(^ ^ are consistent with the experimental data in overdoped region. 



V. SUMMARY AND CONCLUDING REMARKS 



In summary, we have studied a model of HTSC proposed by C. M. Varma which predicts QCP associated with the 
circulating current order parameter. The results derived with the ladder approximation indicates that the system is like 
Fermi liquid near the quantum critical point and in the disordered regime. In particular, the imaginary part of electron 
self-energy is proportional to [max(e,T)]^ up to minor logarithmic correction. Thus, the ladder approximation is not 
sufficient in describing the properties anomalous normal state, for which MFL theory gives adequate phenomenological 
description. For the proper explanation of the anomalous normal state properties the consideration of higher order 
corrections beyond ladder approximation seems to be crucial. 

To get some clues on the role of higer order correction, let us estimate a higher order digram which corresponds 
to the vertex correction of the susceptibility D22(*^,q) beyond ladder approximation. Let us consider the renor- 
malization of the susceptibility by the 4th order terms of ([III). The full momentum and frequency dependences 
of I{c{(li)^^D{<li)^^E{<li) should be retained in evaluating the higher order correction. As has been discussed 
in Sec. IV, the coupling between the anti-bonding conduction electrons and the order parameter is strongly sup- 
pressed for small momentum transfer. As can be shown directly by substituting (A6) into (^) and (^), the coupling 
^p+q^p-^2(q) + H.C. is not suppressed at small q. Thus, we have to consider a diagram which renormalize D22 with 
only (ap^q6pA'2(q) -l- II.C.)-type vertex. The first of such diagram is (See Fig.l), 
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FIG. 1. A higher order diagram which renormalize the susceptibility {N2N2). The roman letter A, B indicates the Green 
function Ga/b (see Eq.(A7).), and the wiggly line is the propagator D22{k). 



With the notations q ~ {iuj, q),p = {ie, p), k = {iil, k), the analytic expression of Fig. 1 is 

(Fig.l) --V^Y. Gb{p)Ga{p + k)GB{p + k + q)GA{p + q)D22{k), (40) 

where the non-singular parts of vertices are omitted. The most singular contribution coming from the low energy 
region can be estimated as follows: First, put fc = in Ga/b and do p integral. Second, carry out the remaining 
integral in fc. If any infrared divergence occurs, cut the divergence by max(q, 62)- The result of p integral is 



/p ~ ^ In 



1 



{EBikF)-^i)^ 



(41) 



Major contribution in momentum integral comes from |p| > pp, where q can be neglected. Usually oj^ <^ 
{EBikp) — , then Ip is non-singular. The diagram Fig.l gives non-singular contribution, which just renormal- 
izes some coefficients of the susceptibility. Therefore, the straightforward perturbative computation at two loop, order 
does not give indications of singular behaviour. However, if we assume, following Varma's heuristic derivationo, that 
in the strong coupling limit (large V) Esikp) — fi can be replaced by an excitonic energy scale which corresponds S2 
in our notation, then Ip becomes highly singular as we approach the quantum critical point (^2 ~^ 0). The justification 
of the above replacement and the further investigation of consequenp^ of the singularity will be further studied. The 
approach using the (non-perturbative) dynamical mean-field theoryfij is currently under investigation. 
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APPENDIX A: DIAGONALIZATION OF KINETIC TERM 



The fermion Green functions can be obtained by inverting the kernel matrix Kq^ Eq.(^. The diagnonalization of 
the free Hamiltonian Hq gives three bands: the anti-bonding i?yi(k), the bonding EbO^), and the non-bonding EcO^)- 
In the limiting case t' — xtpp = 0, the closed expression of each band is 



EAA^)-"4±^,Ec{k)^0, 



(Al) 



1/2 



where -B(k) = (cd)^ + (4t)^s^(k) with s^(k) = s^(k) + s^(k). Finite t' gives dispersion to the non-bonding 
band which we ignore. However, t' plays a crucial role in the description of the current pattern of the ordered phase. 
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Since only the disorded phase is considered in this paper, we will take t' = 0. Let U be the unitary matrix which 
diagonahzes Kq^. 



■ Ea{^) - H-ie 

i;s(k)-/i-ie 

Ec{^)-ti-ie, 



(A2) 



Kq ^ . The matrix of fermion Green's function is 



The explicit form of the unitary matrix is given hy U = (|ri), |r2), Irs)), where |rj) is the j-th eigen-column vector of 

V 



G = U 



ie+/i-_EA(k) 







1 

ie+(Li-BB(k) 





(A3) 



ie+([t-Bc(k) 



where the row and the column of G is indexed by d,Px,Py The explicit form of the unitary matrix U 



U 



I gA(k) gB(k) n. \ 

I a(k) $B{k) U ^ 

2tSa; (k) 2ta^ (k) Sy (k) 

lUkF CB(k) fclkj 

. 2ts^(k) 2fa^(k) -s^(k) , 

\ a(k) «B(k) |c(k) / 



(A4) 



where U/b{^) = E^^j^ik) + {2tfs\k) 



1/2 



u 



When Cd = 0, the matrix elements of U simplify considerably. 



/ J_ L 



V2s V2s 
\ V2s V2s 



(A5) 



The relation between the diagonalized fermion operators a,b,c and the original {d,Px,Py) operators are 

[a,6,c]*■^ = ^/t[d,p„p,]*^ (A6) 
where tr denotes the tranpose of a matrix. In case = 0, the direct relations between fermion Green functions are 



1 ^ 

Gdd = 2^^^^ ^b), Gdp^ = ^(^-4 - Gb), Gdpy = '^i^A - Gb), 



2 g2 g2 2 

^PxP^ = ^(*^^ + + -J-Gc, Gp^py = 7^{Ga + Gb) + -J-Gc, Gp^py = -^{Ga + Gb - 2Gc), 



2s2 



(A7) 



where Ga/b/c = ,e+M-i;l/B/c(k) 



APPENDIX B: POLARIZATION FUNCTIONS 



The explicit expressions of the polarization functions llaiq) are 



^A{q) = -^J2Y1 '^ddip + q) [slip) Gp^pM + 4(P) Gpypyip) - 2s,(p)sj,(p) 

Gp^Py (P) 

" e p 

IlB{q) = -j^'^Yl [GdpAp)sx{p) - Gdpy{p)sy{p)j \Gdp^{p + q)sx{p + (i) - Gdp^ip + q)sy{p + (i) 

nc(gi) = ^ GddiPl) GddiPs) (4(P2)Gp,p,(p2) + 4(P2)<^P,/P«(P2) - 2Sx(p2)%(P2)Gp^p^(P2)), 

P 

X (s^(p4)Gp,p, (P4) + 4 (P4)Gp„p„ (P4) - 2s^(p4)Sj,(p4) Gp^p^ (^4)) , 
nD(gi) = 21/" ^ Gdd(pi) (sl{p2)Gp^p^{p2) + sl{p2)Gp^py{p2) - 2s^(p2)Sj;(p2)Gp^p^(p2)) 
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GdpAP3)SxiP3) - Gdpy{P3)Sy{p3) GdpAp4)SxiPi) - Gdpy{pi)Sy{p4:) 



^E{qz) = ^ n [GdpAPt)Sx{Pt) - Gdpy{p^)Sy{p^ 



1=1,4 



(Bl) 



where the signs of momentum and energy along the loop of the 4th order d iag rams are to be implicitly understood. 
Nq is the number of lattice sites. Let us assume = 0. Using the relation (A7) we can compute IlA{q) and IlB{q)- 



(4(p)-4(p))^ 

2s2(p) 

, (4(p)-4(p))^ 



(4(p)-4(p))' 

2s2(p) 



■GA{p + q)GA{p) 



{Ga{p + q)GB{p) + Ga{p)Gb{p + q)) 



2s2(p) 



GB{p + q)GB{jp) 



44(p)4(p) 

S2(p) 



GB{p + q)Gc{p) 



4s2(p)s2(p) 



(B2) 



''^''^ i (27r)2 A. 2s2(p) 2s2(p + q) 

X \Ga{p)Ga{p + q) + Gb{p)Gb{p + q) - Ga{p)Gb{p + q) - Gb{p)Ga{p + q) 



(B3) 



The angular factors do not introduce any singular features in the integral, and they will be ignored. The calculation 
of GaGa correlation function gives the well-known result (at zero temperature), 



d^p 

{2nf 



tY,Ga{p)Ga{p + 9) ^ ~N{Ef) 



1 



vfM 



(B4) 



where N{Ep) is the unrenormalized density of states at the Fermi level, and the purely numerical constants are 
suppressed. The result (B4) is valid for |q| < pp, \oj\ < w_F|q|. In the regime, |w| ^ w_F|q|, the correlation function 
( p^ is proprotional to [qP/w^^ therefore, it is small. The simrmation over frequency of GAip + q)GB{p) correlation 
function gives 



(Pp npiEBip) - /i) - np{EAip + q) - /i) 



iuj + EB{p)-EAip + q) 



Similary, 



Pp np{EAip) - /i) - np{EBip + q) - m) 



-iLU + Ea{p) - Eb{p + q) 



(B5) 



(B6) 



At low temperature np{EB{p) — /i) ~ 1, therefore, the anti-bonding electrons should lie outside the Fermi surface. 
Thus, for the momentum |q| < pp, ( p5| ) and (B6) do not depend on q. Carrying out the integrals, we get for < 2t, 



/j2 
J^TJ2 {GBip)GAip + q) + Ga{p)Gb{p + q) 



Aq - Pf 



{2ty 



Ao is the momentum cut-off, or the size of Brillouin zone. In the opposite limit, |a;| S> 2t, 



/j2 
(^^E [Gb{p)Ga{p + q)+ Ga{p)Gb{p + q) 



27r \uj\ 



GaGc can be computed analogously. 

_d 

(2 



—T2_^Gc{p)GA{p + q) z In- 



2tt 



(B7) 



(B8) 



(B9) 
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where is the energy cut-off. In the Hmit, < /i, 

— N{Ef) 



(BIO) 



The correlation functions GbGb and GbGc give negligible contributions at low temperature because the bonding 
band and the non-bonding band are fully occupied. Combining all of the above results, we can write 



Ha (9) = hN{Ep)V^ 
ILB{q)=b2N{EF)V^ 



UJ 



Cl 



C2 + |qr + 



OJ 



vF\q\ m 

2 



In 



— lUJ 



Co +M 



CO 



(Bll) 



where bi,b2,ci are the positive numerical constants. C2 is also a numerical constants which can positive or negative 
depending on doping and bandwidth. In case of our interest, C2 is positive. 
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